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Kubo's formula relates bulk viscosity to the retarded Green's function of the trace 
of the energy-momentum tensor. Using low energy theorems of QCD for the latter we 
derive the formula which relates the bulk viscosity to the energy density and pressure 
of hot matter. We then employ the available lattice QCD data to extract the bulk 
viscosity as a function of temperature. We find that close to the deconfinement 
temperature bulk viscosity becomes large, with viscosity-to-entropy ratio C/s ~ 1- 

One of the most striking results coming from RHIC heavy ion program is the observation 
that hot QCD matter created in Au — Au coUisions behaves hke an almost ideal liquid 
rather than a gas of quarks and gluons . Indeed, hydrodynamical simulations ol' 

nuclear collisions at RHIC (see e.g. 0, 7]) indicate that the shear viscosity of QCD plasma 
is very low even though a quantitative determination is significantly affected by the initial 
conditions j^. This observation does not yet have any theoretical explanation due to an 
enormous complexity of QCD in the regime of strong coupling. This is why the information 
inferred from the studies of gauge theories treatable at strong coupling such as = 4 SUSY 
Yang-Mills theory is both timely and valuable. The study of shear viscosity in this theory 
using the holographic AdS/CFT correspondence has indicated that the shear viscosity t] at 
strong coupling is small, with the viscosity-to-entropy ratio not far from the conjectured 



bound of rj/s = I/Att 



However = 4 SUSY Yang-Mills theory is quite different from QCD; in particular it 
possesses exact conformal invar iance whereas the breaking of conformal invar iance in QCD is 
responsible for the salient features of hadronic world including the asymptotic freedom Uj , 
confinement, and deconfinement phase transition at high temperature^. Mathematically, 
conformal invariance implies the conservation of dilatational current s^: = 0. Since 

the divergence of dilatational current in field theory is equal to the trace of the energy- 
momentum tensor (9^s^ = 6^, in conformally invariant theories 6^ = 0. In QCD, in the 
chiral limit of massless quarks the trace of the energy-momentum tensor is also equal to zero 
at the classical level. However quantum effects break conformal invariance isl. [l^: 

= 0^, = J2 + f^Tr G^^G,, , (1) 

where (3{g) is the QCD /^-function, which governs the behavior of the running coupling: 

= /3(s); (2) 

note that we have included the coupling g in the definition of the gluon fields and have not 
written down explicitly the anomalous dimension correction to the quark mass term. 

How would this breaking of conformal invariance manifest itself in the transport properties 
of QCD plasma? How big are the effects arising from it? The transport coefficient of the 
plasma which is directly related to its conformal properties is the bulk viscosity; indeed, it is 
related by Kubo's formula to the correlation function of the trace of the energy-momentum 
tensor: 

C = - hm - / dt d're^-' mx),e^^iO)]) . (3) 

It is clear from ([3]) that for any conformally invariant theory with 6'^ = 6* = the bulk 
viscosity should vanish. 

The perturbative evaluation of the bulk viscosity ( of QCD plasma has been performed 
recently [l^, and yielded a very small value, with C/-^ ~ 10^'^. The parametric smallness 
of bulk viscosity can be easily understood from eqs ([3]) and ([T]) which show that ( ^ a^, in 



The effects of conformal symmetry breaking on bulk viscosity of SUSY Yang-Mills theory have been 
studied in the framework of the AdS / CFT correspondence in ref (l^ . 



accord with the result of ref. [15|. This would seem to suggest that bulk viscosity effects in 
the quark-gluon plasma are unimportant. However, perturbative expansions at temperatures 
close to the critical one are not applicable, so at moderate temperatures one has to rely on 
lattice QCD calculations. Lattice calculations of the equation of state become increasingly 



precise; however, the direct calculations of transport coefficients 



lave been notoriously diffi- 



cult. Two calculations have been reported for shear viscosity |16l. Il7l|. including a recent high 
statistics study [17]. Both indicate that rj/ s is not much higher than the conjectured bound 
of l/47r; no lattice calculations of the bulk viscosity have been reported so far. Fortunately, 
the correlation function of the trace of the energy-momentum tensor in QCD is constrained 
by the low-energy theorems, which do not rely on perturbation theory. They can thus be 
used to express the bulk viscosity in terms of the "interaction measure" (6) = S — 3P where 
S is the energy density and P is the pressure, which are measured on the lattice with high 
precision. Such a study is the subject of this Letter. 

The calculation of the bulk viscosity starts with the Kubo's formula (jSj) (we follow the 
definitions and notations of [28]). Introducing the retarded Green's function we can re- write 
(El) as 

C = - lim - / dt d^r e*^* iG^(x) = - lim - iG^(cu, 0) = — lim - ImG^ico, 0) . (4) 
9 uj-^o CO Jq J 9 ^~*o uj 9 i^^o CO 

The last equation follows from the fact that due to P-invariance, function ImG^(a;, 0) is odd 
in UJ while ReG^(ti;, 0) is even in uj. Let us define the spectral density 

p{uj,p} = --lmG^{u;,p}. (5) 

TT 

Using the Kramers- Kronig relation the retarded Green's function can be represented as 



G'^{uj,p) = - ^-^du= '^\ du (6) 

^ J-oD u — Lu — ie J„oo UJ — u + ie 

The retarded Green's function G^{uj,p) of a bosonic excitation is related to the Euclidean 
Green's function G^{uj,p) by analytic continuation 

G^{u,p) = ~G^{iuo,p) , u; > 0. (7) 

Using (j6]) and the fact that p{uj,p) = —p{—uj,p) we recover 

G^(0,0) = 2 r P^^du. (8) 
Jo u 
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As we discussed above, the scale symmetry of QCD lagrangian is broken by quantum 
vacuum fluctuations. As a result the trace of the energy momentum tensor 9 acquires a 
non-zero vacuum expectation value. The correlation functions constructed out of opera- 
tors 9{x) satisfy a chain of low energy theorems (LET) which are a consequence of the 
renormalization group invariance of observable quantities [3]. These low-energy theorems 
entirely determine the dynamics of the effective low-energy theory. This effective theory has 
an elegant geometrical interpretation jl9]; in particular, gluodynamics can be represented 
as a classical theory formulated on a curved (conformally flat) space-time background {20] . 
At finite temperature, the breaking of scale invariance by quantum fluctuations results in 
9 = £ — 3P 7^ clearly observed on the lattice for SU (3) gluodynamics 2l| : the presence of 



large Nc 



quarks 22l| i ncluding the physical case of two light and a strange quark [23], or considering 



24| does not change this conclusion. 



The LET of Ref. [iSl, \l^ were generalized to the case of finite temperature in (26|, l2 
The lowest in the chain of relations reads (at zero baryon chemical potential): 

G^(0, 0) = 1 d'x {T9{x), 9(0)) = (t-^ - 4^ {9)t • (9) 

To relate the thermal expectation value of {9)t to the quantity {S — 3P)lat computed on 
the lattice, we should keep in mind that 

{S - 3P)lat = {0)t - {0)o, (10) 

i.e. the zero-temperature expectation value of the trace of the energy-momentum tensor 

{9)o = -4|e,| (11) 

has to be subtracted; it is related to the vacuum energy density < 0. Now, using ([8]), ([9]) 
and (fTOj) we derive the following sum rule 

Jo u V or) ^ dT T4 ^ I ''I ' y ) 

This exact relation is the main result of our paper. 



In order to extract the bulk viscosity C, from ( 1121) we need to make an ansatz for the 
spectral density p. At high frequency, the spectral density should be described by perturba- 
tion theory; however, the perturbative (divergent) contribution has been subtracted in the 
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FIG. 1: The ratio of bulk viscosity to the entropy density for SU{3) gluodynamics. We have used 
|e^,| = 0.62 and = 0.28 GeV [21]. 



definition of tlie quantities on tlie r.fi.s. of the sum rule f[T2|) . and so we should not include 
the perturbative continuum^ p{u) ~ on the l.h.s. as well. In the small frequency region, 
we will assume the following ansatz 

which satisfies ^ and (jS])- Substituting f[T5]) in f[T^ we arrive at 

S(£-3PW^ieMl. (14) 

A peculiar feature of this result is that the bulk viscosity is linear in the difference £ — 3P, 
rather than quadratic as naively implied by the Kubo's formula. This is similar to the 
strong coupling result obtained for the non-conformal supersymmetric Yang-Mills gauge 



plasma [l2|] 



The parameter uq = ujq{T) is a scale at which the perturbation theory becomes valid. 
On dimensional grounds, we expect it to be proportional to the temperature, o^i ~ T. We 



estimate it as the scale at which the lattice calculations of the running coupling 29|] coincide 
with the perturbative expression at a given temperature. In the region 1 < T/Tc < 3 we 
find LJ ^ (T/Tc) 1.4 GeV. Now we are ready to use (|T^ to extract the bulk viscosity from 
the lattice data. 



^ For an explicit perturbative expression and a discussion of the properties of p{u) at small frequencies see 
e.g. H- 
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The results of the numerical calculation using as an input the high precision lattice data 



2l| are displayed in Fig. [H One can see that away from the bulk viscosity is small, in 



accord with the expectations based on the perturbative results 15|]. However, close to Tc the 
rapid growth oi S — 3P causes a dramatic increase of bulk viscosity. Basing on the lattice 



results 



16 



17( 1 which indicate that the shear viscosity remains small close to Tc, we expect 



that bulk viscosity will be the dominant correction to the ideal hydrodynamical behavior in 
the vicinity of the deconfinement phase transition. 
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